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AN EXACT CLOSED-FORM SOLUTION FOR CONSTANT-AREA COMPRESSIBLE 
FLOW WITH FRICTION AND HEAT TRANSFER 
by Jonas I. Sturas 
Lewis Research Center 

SUMMARY 

The well-known differential equation for the one-dimensional flow of a compressible 
fluid with heat transfer and wall friction has no known solution in closed form for the 
general case. This report presents a closed form solution for the special case of con- 
stant heat flux per unit length and constant specific heat. The solution was obtained by 
choosing the square of a dimensionless flow parameter as one of the independent vari- 
ables to describe the flow. From this exact solution, an approximate simplified form is 
derived that is applicable for predicting low Mach number flow performance character- 
istics for many types of constant area passages in internal flow. The data included in 
this report are considered sufficiently accurate for use as a guide in analyzing and de- 
signing internal gas flow systems. 


INTRODUCTION 

In this report one-dimensional flow characteristics of a compressible fluid under 
simultaneous action of friction and heat transfer in a constant area passage are investi- 
gated. The pressure drop in a constant area passage results from the simultaneous ac- 
tion of friction and heat transfer. In the technical literature, the various equations 
describing the flow characteristics of a compressible fluid under the simultaneous influ- 
ence of friction and arbitrary heat transfer are not expressed in closed form (see, e.g. , 
section 8.9 of ref. 1). Presently, no generalized integrating factor has been found that 
would permit closed-form integration for any fluid temperature variation. 

In the present analysis, a completely integrated solution is presented for the case 
in which the total temperature of the fluid varies linearly with distance through a con- 
stant flow area passage. For constant specific heat, this solution constitutes flow at 
constant heat flux. From this exact solution, an approximate simplified form is derived 



that is applicable for predicting one -dimensional flow performance for most types of 
operation in constant area passages in one -dimensional flow. 

This analysis develops a working chart and relations that simplify the determination 
of the Mach number change in cooling and heating passages and, consequently, simplifies 
the calculation of the pressure drop required across these passages to obtain the neces- 
sary mass flow of gas. The method developed applies to the flow of gas through straight 
ducts having any constant shape, cross-sectional area, and roughness. 

The solution is presented in the form of working charts constructed for a ratio of 
specific heats equal to 1. 40 and 1. 30. The variables cover ranges sufficiently large to 
include many applications of engineering interest. The data included in this report are 
considered to be sufficiently accurate for use as a guide in analyzing and designing inter- 
nal gas flow systems for aircraft. 


SYMBOLS 

A duct flow area, m^; ft^ 
a speed of sound, m/sec; ft/sec 

c specific heat at constant pressure, J/kg-K; Btu/lbm-°R 

r 

c y specific heat at constant volume, J/kg-K; Btu/(lb/°R) 

D equivalent hydraulic diameter, m; ft 
Fj friction force, N; lbf 

o 

f local friction factor = r/pv /2, dimensionless 

W 

2 2 
g c conversion factor, (kg-m)/(N-sec ); (lbm-ft)/(lbf-sec ) 

J Joule's constant, N-m/J; lbf-ft/Btu 

k heat flux parameter, q/CpT q (see eq. (2)), dimensionless 
M Mach number in duct at distance x, dimensionless 
P total pressure at distance x, N/m ; lbf/ft 
p static pressure at distance x, N/m ; lbf/ft 
Q heat input to stream, J/kg; Btu/lbm 
q constant heat flux, J/kg; Btu/lbm 

R gas constant, N-m/kg-K; (lbf ft)/(lbm-°R) 

T total temperature at distance x, K (°R) 
t static temperature at distance x, K (°R) 
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I 


V dimensionless flow parameter at distance x, (v/ < y/g c RT 0 ), dimensionless 
v axial velocity in flow passage distance x, m/sec; ft/sec 
w weight flow rate of fluid, kg/sec; lbm/sec 
x axial distance through flow passage, m; ft 
y ratio of specific heats, Cp/c y , dimensionless 
A increment of a dimensionless parameter 

(i integrating factor, dimensionless 


s 


friction-distance parameter 



o o 

p mass density, kg/m ; lbm/ft 

r defined by eq. (2) as ratio (T/Tq), dimensionless 

O O 

r w frictional shearing stress at the wall, N/m ; lbf/ft 
Subscripts: 

D duct 

l limiting condition 

0 at x = 0 


METHOD OF ANALYSIS 

For one-dimensional flow, the pressure at any point in a fluid is uniquely defined by 
the mass velocity (G = w/A D = pv), the temperature, and the Mach number at that point. 
The following development deals with the variation of Mach number along a constant area 
duct. If the variation of Mach number is known, the corresponding pressure distribution 
along the duct is readily determined from the continuity equation for the one -dimensional 
flow of a perfect fluid. 

The investigation of frictional effects will be confined to steady one-dimensional flow 
of a compressible fluid in a constant cross-sectional area duct. For these conditions a 
flow model and the forces acting on a fluid element in the x-direction are illustrated in 
the following sketch: 
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6F, 


Mass of fluid element 

/- <f m ■ *p'g c )A D dx 


f / 


v <r 


(pA D )- 

V“ 


i r-n 


i i 


-(p + dp)A D 
- (v + dv) 


-x dx— 


(a) 


or 


F x = pA D “ (p + d P) A D " 6F f (!) 

The mechanisms of friction and heat transfer are so similar that one cannot be had 
without the other (ref. 1, p. 242). Therefore, in order to calculate the pressure drop 
and other flow characteristics more accurately, it is necessary to take account of fric- 
tion as well as of heat transfer. In this analysis we shall assume the flow area, specific 
heat, the weight rate flow, and the heat flux all to be constant. Therefore, the rate of 
heat transfer may be expressed by means of the energy equation in terms of the increase 
in stagnation enthalpy. Thus, 


dQ - c p dT 
= c dT 

\dU P 

M = T - T 0 



q 

c p T 0 


I 


or 


r = 1 + k£ 


(2) 
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f 


where 


T 



and 


k = ^L- 

c p T 0 

Therefore, the value of k is primarily determined by the heat flux q. Also, from 
equation (2) it can be seen that the assumptions of constant specific heat and constant 
heat flux are directly analogous to a linear variation in stagnation temperature T along 
the duct. 

If the dimensionless parameters V, and t/Tg are defined by the following ex- 
pressions 


V s 


v 

Vg c RT o 


(3) 


£ 3 — /' X f dx = (4) 

D S Q D 

_L = T_ _ y - l/ v 2 \ 

T 0 T 0 Uc RT 0/ 

= r - V 2 (5) 

2y 

where j = T/Tq, then the relation between the standard definition of Mach number 
(M = v/a) and the dimensionless flow parameter V results: 


M 


2 



( 6 ) 
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or 


v 2 = 


r M 2 T 


1 + 


A) 


( 7 ) 


M 


For assumed conditions, the one -dimensional steady-state momentum equation for a 
compressible flow given in reference 1 (eq. (8. 21)) can be rearranged to this form 


dv 2 + 2g c R^dt + t = 0 


Making proper substitutions for friction-distance parameters £ and axial velocity v, 
the momentum equation becomes 

dV 2 + 2 (— + — • + V 2 d| = 0 

\ T 0 T 0 p) 

Logarithmic differentiation of the continuity equation for constant flow area gives 

dp _ dv _ dV 

p v V 

Eliminating the density term results in a new expression of the momentum equation: 


dv 2 + 2[— - — — ) + V 2 d£ = 0 
\ T 0 T 0 V / 


(8) 


A combination of equations (5) and (8) yields 


V 2 + 2 — )d£ = ( — — 


y±l)dv : 

2y 


(9) 


Substitution of equation (2) into equation (9) results in 


d* - 




V 2 (V 2 + 2k) 


dV 2 = 


_ y+A v 2 \ 

2y | dV z 
V 2 + 2k / V 2 


( 10 ) 
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Equation (10) is linear in and, by multiplying through with the integrating factor, 


M = 


< V_ + 2k ' 
( V 2 ) 


1/2 


which satisfies the necessary requirements that 


— (juM) si-(juN) 


av" 


31 


(ID 


where 


M = 1 


-(k^+l-^V 2 ) 

n = — - ^ — y 

V 2 (V 2 + 2k) 

(see ref. 2), equation (10) becomes an exact differential equation, which in mathematical 
notation can be written as 


f(V 2 )d| + £g(V 2 )dV 2 + h(V 2 )dV 2 = 0 


(12) 


Integration of equation (12) gives 



y + 1 

y 


2 V 2 

(V^ + 2k) 


In 


Vq + k + V Q (v 2 + 2k 


1/2 


2 2 1/2 

V" 1 + k + V(V" ! + 2k) 


(13) 


In order to solve for % for adiabatic conditions, the following steps were taken to 
eliminate the heat flux parameter k from equation (13): 
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Substitution of equation (7) into equation (14) with k = 0, gives 





( 15 ) 
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It is interesting to note that equation (15) is of the same form as a combination of 
equations (6. 20) and (6. 21) given on page 167 of reference 1 for adiabatic flow with fric- 
tion. 

Useful relations of dimensionless weight rate flow parameters can be written as 
follows: 

(1) In terms of static pressure: Substitution for M from equation (6) into the fol- 
lowing equation taken from reference 1 (p. 82), 



which expresses the weight rate flow parameters in the classical form, results in a 
dimensionless weight rate flow parameter for diabatic conditions 


V [j / y- 1\ vfl 
A D g c p T 1 / 2 L \ 2y / T J 


(16) 


(2) In terms of total pressure: Again, substitution for M from equation (6) into 
the following equation, taken from reference 1 (p. 84) 


wVt 


a d p 



(y+l)/[2(y-l)] 


M 


results in a dimensionless weight rate flow parameters for diabatic conditions 


wVe c RT 

a dSc p 



(17) 


From equations (16) and (17), it follows that the static and total pressure ratios for 
diabatic conditions will be, respectively, 


JP_ 

P0 



'1 + 




>^M 2 


1(1 + k|) 


1 + 


1/2 


(18) 
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and 


P. = /^\l + k| )l/2| 
P 0 \ M / 






,1 + 




Dividing equation (18) or (19), gives 




^v-rAr-i) 

M o/ 


( 19 ) 


( 20 ) 


Equation (20) indicates for subsonic flow that total pressure drop always exceeds 
static pressure drop in a constant area compressible flow because of fluid acceleration. 

An approximate explicit expression for V may be obtained by use of a Maclaurin 
expansion of the logarithmic quantity that appears in equation (13). Using the first two 
terms of such an expansion eventually leads to the result 


V = (1 + k|) 




0.0001 k|) 1/2 |l - V 2 | j| 


y + J 

12y 


-j ^20 Vq + 15k^ + 2 + k£ j 


- 1/2 


( 21 ) 


LIMITING CONDITIONS FOR DIABATIC FLOW 

In compressible subsonic flow, as the gas flows along the duct, its static pressure 
decreases while its velocity increases. As shown in reference 1 (pp. 255 to 256), the 
effect of the wall friction, for diabatic (nonadiabatic ) conditions, is to accelerate the 
flow velocity. 

Equation (9) can be rearranged to the following form 



k£ + 1 - | 2 

V 2 (V 2 + 2k) 


When d£ = 0, there can be no change in Mach number; that is, 
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k£ + 1 
d| 

dV 2 V 2 (V 2 + 2k) 


- (y+ 1 

V 2y , 


V 


= 0 


from which, the limiting parameter is 


v ' - u + 


r2 ; 


Substituting into equation (6), gives 


Mf = 


V 


y(l + k^) - vf 


2 y - V 




'y-i 


= 1. o 


( 22 ) 


Therefore, 


=1.0 (23) 

It is interesting to note that, for a constant heat flux, as is the case assumed in this 
analysis, the same result (that is, = 1. 0) could be obtained using equation (8. 71a) of 
reference 1. 

The relation between the limiting or choking static pressure ratio P^/p Q and the 
friction-distance parameter £ can be obtained by substituting =1.0 into equa- 
tion (18). Thus, 



(24) 
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and for adiabatic conditions when k = 0 equation (24) reduces to equation (6. 22) of ref- 
erence 1. Also, the relation between the limiting or choking total pressure ratio P^/Pq 
and the friction distance parameter 4 can be obtained by substituting M = 1. 0 into 
equation (19): 



= (M 0 )(l + k£) 1/2 



(y+l)/[2(y-l)] 


(25) 


Again, for adiabatic conditions equation (25) reduces to equation (6. 26) of reference 1. 

As can be seen from equations (25) and (24), the heat addition to the stream tends 
to increase choking pressure ratio P^/Pq and P^/P q» respectively, by the factor of 
d + kl) 1 / 2 . 


RESULTS AND DISCUSSION 
Diabatic Flow with Friction 

The mechanisms of friction and heat transfer are similar because friction and heat 
transfer cannot exist independently of each other. Therefore, in order to calculate the 
pressure drop and other flow parameters of interest, for example, such as bled hot or 
cold gas flow from a turbofan engine, it is necessary to consider both friction and heat 
transfer together. 

In this analysis the heat transfer is accounted for by assuming a constant heat 
flux q. As has been proved under the method of analysis (see eq. (2) and fig. 1), a con- 
stant heat flux is directly analogous to the case of a linear variation of total temperature 
along the duct. 

For one -dimensional flow, the pressure at any point in a fluid is uniquely defined 
by the mass velocity (w/A^ = pv), the temperature T, and the Mach number M at that 
point. Once the local Mach number along the duct is known, the other flow parameters 
such as total or static pressure drop, weight rate flow per unit area, and the tempera- 
ture distribution along the duct are readily determined. Figure 2 represents the varia- 
tion of local Mach number M against the friction-distance parameter £ for an inlet 
Mach number of 0. 15, y - 1.3 and 1. 4, and a range of heat flux parameters 
0 < k < 0. 3. 

As can be seen from this plot, for a given friction-distance parameter £, the 
Mach number is very sensitive to an increase in heat flux parameter k. 
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0 4 8 12 16 20 24 28 

Friction-distance parameter, £ 


Figure 1. - Linear variation of total temperature with friction- 
distance parameter. 


, ■ -x- „ u Heat flux 

Limiting Mach parameter, 



Friction-distance parameter, 


Figure 2. - Variation of local Mach number in the duct with friction- 
distance parameter. inlet Mach number, 0.15. 
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Friction -distance parameter, 


Figure 3. - Static pressure ratio as function of friction-distance param- 
eter. Local Mach number, 0. 15; ratio of specific heats, 1.4. Limit- 
ing condition is determined by choked Mach number for diabatic con- 
ditions. 



Friction-distance parameter, £ 


Figure 4. - Total pressure ratio as function of friction-distance param- 
eter. Local Mach number, 0. 15; ratio of specific heats, 1.4. Limiting 
condition is determined by choked Mach number for diabatic conditions. 


Figure 3 presents the ratio of the local to inlet static pressures as a function of 
friction-distance parameter £, and figure 4 presents the ratio of the local to inlet total 
pressure as a function of the friction-distance parameter In figure 5, the weight 
flow parameter w^T/A^P is plotted against friction-distance parameter ij. Figures 
3 to 5 are presented for a range of heat flux parameters 0 < k < 0. 4 that were calcu- 
lated using equation (2) and from the test data presented in reference 3. In all of the 
figures, the inlet Mach number equals 0. 15 and y = 1. 4. The dashed line in the figures 
represents the limiting or choking conditions, defined by equations (24) and (25), re- 
spectively. 

If the variation of local Mach number is known, the corresponding corrected weight 
flow rate for diabatic conditions is readily determined from equation (16) in terms of 
local static pressure or from equation (17) in terms of the local total pressure along 
the duct. 


Heat flux 
parameter, 
k 



Figure 5. - Variation of weight flow parameter with friction distance. 
Inlet Mach number, 0. 15; ratio of specific heats, 1.4. 
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Exact local duct Mach number, M 


Figure 6. - Percent error introduced in local Mach number by using 
approximate solution (eq. (2D). Ratio of specific heat, 1.4. 

Accuracy of Approximate Analysis 

Figure 6 presents the relative error of the local Mach number calculated using an 
approximate equation (21) and an exact equation (13). It can be seen that, for an inlet 
Mach number equal to 0. 1, an approximate equation (21) is accurate to within ±1. 5 per- 
cent relative to an exact equation (13) up to a local Mach number of 0. 3 for all values of 
heat flux parameters k. However, for an inlet Mach number of 0. 2, an accuracy of 
approximately ±2 percent is maintained up to a local Mach number of 0.5. Although it 
is not shown, the error in the local Mach number is insensitive to a decrease in y from 
1. 4 to 1. 3. 

In general, it can be seen from figure 6 that, for all values of heat flux parame- 
ters k, the relative error of the local Mach number calculated using an approximate 
equation (21) and an exact equation (13) progressively increases with an increase in inlet 
Mach number Mq. 

SUMMARY OF RESULTS 

For constant area flow, in which the heat flux and specific heat remain constant, a 
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method is developed for calculating the local flow Mach number, the pressure drop, and, 
therefore, the weight flow rate of a compressible fluid under the simultaneous action of 
friction and heat transfer. The results are presented in the form of relations and work- 
ing charts. 

A closed form solution to the problem is achieved by choosing the square of a di- 
mensionless flow parameter as one of the independent variables to describe the flow. 
From this exact solution, an approximate simplified solution is presented which is in 
good agreement with the exact solution for low subsonic Mach numbers. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 21, 1971, 

764-72. 
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